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Abstract. Recently, angle-resolved photocmission spectroscopy (ARPES) has 
been used to highlight an anomalously large band renormalization at high binding 
energies in cuprate superconductors: the high energy "waterfall" or high energy 
*vj , anomaly (HEA). This paper demonstrates, using a combination of new ARPES 

measurements and quantum Monte Carlo simulations, that the HEA is not simply 
the byproduct of matrix element effects, but rather represents a cross-over from a 
quasi-particle band at low binding energies near the Fermi level to valence bands 
at higher binding energy, assumed to be of strong oxygen character, in both 
hole- and electron-doped cuprates. While photocmission matrix elements clearly 
play a role in changing the aesthetic appearance of the band dispersion, i.e. the 
^^, ■ "waterfall" -like behavior, they provide an inadequate description for the physics 

f^ ' that underlies the strong band renormalization giving rise to the HEA. Model 

QQ ' calculations of the single-band Hubbard Hamiltonian showcase the role played 

^^v ' by correlations in the formation of the HEA and uncover significant differences 

in the HEA energy scale for hole- and electron-doped cuprates. In addition, 
this approach properly captures the transfer of spectral weight accompanying 
- both hole and electron doping in a correlated material and provides a unifying 

J^ ' description of the HEA across both sides of the cuprate phase diagram. 

.^, 

PACS numbers: 79.60.-i, 71.10.Fd, 74.25.Jb, 74.72.-h 
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1. Introduction 

Over the last decade significant advancements have improved the momentum and 
energy resolution of angle-resolved photoemission spectroscopy (ARPES) pT] that 
provides access to the single-particle spectral function yl(k, w). While a significant 
probe of electronic structure in general, such advancements in ARPES have had 
a profound impact on the study of strongly correlated materials. A series of high 
resolution ARPES experiments at binding energies up to ~ 1 — 1.5 eV, made possible 
by these advancements, has revealed the presence of a high energy anomaly (HEA) in 
the dispersion of both hole- [1 H HH [H [3 H] and electron-doped (Sm2_:rCe:rCu04 0, 
Nd2-a;Ce2,Cu04 [TUj, Pri_:rLaCe:rCu04 [TT) high- Tc superconductors. In hole-doped 
compounds the anomaly, at an energy ^ 300 meV, is present in different materials, 
doping values, and phases and more recently has been resolved in electron-doped 
compounds at approximately twice the energy scale. A link for the HEA between hole- 
and electron-doped cuprates can be found in the similar "waterfall" feature observed 
in the half-filled parent insulators [T^] . A number of theories have been advanced to 
explain the presence of this anomaly including spin-charge separation [2] , spin polarons 
[13], in-gap band-tails [Mj, coupling to phonons [4], plasmons [15], or paramagnons 
[3 [m [171 HH] J photoemission matrix elements [7] , and strong correlation or "Mott" 
physics [31 [191 (ZHl [211 [22] ■ These theories include assertions that no HEA should 
be |21| or is [23j present in electron-doped cuprates contrary to other experimental 
evidence [^[iniin]. 

The appearance of the HEA as a "waterfall"- or kink- like feature provides 
conceptual appeal for certain theoretical explanations [H \5\ [151 [Ml [El HH]. However, 
such scenarios based on weak coupling to high energy bosonic modes have severe 
problems. Attributing the HEA to electron-phonon coupling, similar to that used to 
describe the "low-energy kink" in the cuprates j2|l[551[lB], appears to be unhkely given 
the relatively large energy scale and its apparent lack of doping dependence. Coupling 
to paramagnons or spin fluctuations generally satisfies the energy for the HEA, based 
on the superexchange value J ~ 100 meV in cuprates; however, such coupling fails to 
properly address the apparent damping of paramagnon modes with doping [27] and 
cannot account for the dichotomy in energy scales between hole- and electron-doped 
materials. 

In addition, the kink-like appearance gives way to "band break-up" at low photon 
energy [8] or a shallow band dispersion with a characteristic "Y" appearance near the 
zone center in the second Brillouin zone (BZ) Jj, both effects primarily due to changes 
in photoemission matrix elements. Photoemission matrix elements complicate the 
analysis of individual spectra through extrinsic effects [7] , but there remain intrinsic 
band renormalization effects creating a shallow, dispersing quasiparticle band and 
coherent /incoherent crossover at the HEA energy scale [HI [10] for which solely extrinsic 
mechanisms are an inappropriate description |18j . 

Here, new results from ARPES on Ndi.83Ceo.i7Cu04 (NCCO) also reveal the 
presence of a HEA. These results demonstrate changes in the appearance of the band 
dispersion with incident photon energy (photoemission matrix elements) while clearly 
revealing an underlying band renormalization not attributable to matrix element 
effects. The energy scale for the HEA, ~ 500 — 750 meV along the nodal direction, is 
approximately twice that found in the hole-doped high- Tc compounds, as revealed in 
other studies [gifTUlfTT]. 

To understand these results together with previous findings in hole-doped and 
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parent cuprates, we present results from quantum Monte Carlo calculations of the 
single-band Hubbard model, across both sides of the phase diagram and at half-filling. 
Our results indicate that the HEA can be connected to doping and subsequent spectral 
weight transfers into the lower or upper Hubbard band of hole- or electron-doped 
cuprates, respectively. These calculations show that hole or electron doping away from 
the parent insulator promotes the formation of a quasi-particle band crossing Ep , the 
precursors of which can be found in the half- filled parent, and that the HEA represents 
a coherent/incoherent cross-over from this band to valence bands at higher binding 
energies. In addition, the results highlight the electron/hole doping asymmetry of the 
HEA energy scale, in agreement with experimental observations. 

2. Experimental Results 
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Figure 1. (a) Intensity plot along a nodal cut in the first BZ of NCCO (x=0.17) 
taken at 16.75 eV incident energy, (b) Energy distribution curves (EDCs) and 
(c) momentum distribution curves (MDCs) of the cut in (a). The direction of 
the MDCs and EDCs are indicated by the black arrows in panel (a) and curves 
corresponding to kp and Ep are highlighted in each panel. The HEA appears as 
a "waterfall" at an energy scale ~ 500 — 600 meV along this cut. 



ARPES data were taken at beamline 5-4 of the Stanford Synchrotron Radiation 
Lightsource with a Scienta R4000 analyzer at photon energies from 15 to 20 eV. 
The energy resolution was ^ 20 meV with an angular resolution ~ 0.3°. NCCO 
single crystals were grown in 4 atni of oxygen using the traveling-solvent floating-zone 
technique, annealed for 10 h in argon at 970 C followed by 20 h in oxygen at 500 C 
[28] and then cleaved and measured at 10 K and pressures better than 3 x 10~^^ torr. 

Figure 1 displays the ARPES data along a nodal cut ((— 7r/a, —it /a) to (tt/o, 7r/a)) 
for NCCO taken at 16.75 eV. The intensity plot of figure 1(a) clearly shows the HEA 
as a "waterfall" in the dispersion at an energy ~ 500 — 600 meV. This behavior is 
mirrored in the energy distribution curves (EDCs) and the momentum distribution 
curves (MDCs) shown in figures 1(b) and (c), respectively. The EDC and MDC 
derived dispersions are superimposed over the intensity plot of figure 1 (a) . This HEA 
demarcates the transition between a quasi-particle band at low binding energy and 
the oxygen valence bands at energies > 1 eV. While the MDCs appear to show a steep 
dispersion that bends back below the HEA, the EDCs show a shallower band with 
a rapid intensity decrease on the order of the HEA energy scale. This difference is 
similar to that found in the half-filled parent insulators [T^] and hole-doped systems 
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Figure 2. (a) Intensity plots along nodal cuts in tlie first BZ of NCCO (x=0.17) 
taken at tlie indicated incident photon energies. Note tfie aestlietic clianges 
induced by variation in matrix elements with photon energy, (b) EDCs taken 
from the cuts in (a) at the same photon energies. While the HEA appears 
as a "waterfall" along these cuts for several photon energies, the EDC derived 
dispersions indicate a shallow, dispersing band with the HEA energy scale 
~ 500 - 750 meV. 



Elil- 

Figure 2(a) shows a series of nodal cuts at different photon energies highhghting 
the effect of changes to photoemission matrix elements. The data demonstrate 
aesthetic changes to features and apparent negative band velocities below the HEA 
energy scale intimately tied to matrix element effects. However, they also reveal an 
underlying band renormalization producing a HEA on the scale ~ 500 — 750 eV. EDCs 
at these photon energies provide clarity to this point. Figure 2(b) displays EDCs of 
the nodal cuts extracted from the respective data presented in figure 2(a). Each of 
these EDCs shows a shallow dispersing band whose intensity decreases approaching 
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Figure 3. Comparison of tiie HEA between iiole- and electron-doped systems. 
The left panel shows the nodal ARPES data for BSCCO-2201 adapted from [J; 
the right panel shows the nodal data for NCCO at 16.75 eV from this study. 
Matrix element profiles obtained from a 2D fitting procedure 1291 are also shown 
in each panel (the blue arrow indicates that the right most axis should be used for 
matrix element values). While the appearance of the HEA is affected by changes 
in the photoelectric matrix elements, they can not account for the strong band 
renormalization. 



the BZ center, similar to the behavior in hole-doped cuprates [2j[8]. EDC derived 
dispersions yield band bottoms (or the energy where the ARPES intensity falls-ofF) 
in accord with the energy scale set by the HEA, not the deeper energy scales usually 
assigned by LDA estimates and MDC derived dispersions [21 [3] . 

A direct comparison of the ARPES data from BizSrzCuOe (BSCCO-2201) [3] 
and NCCO is shown in figure 3. The HEAs are qualitatively similar in the two 
materials, but there are significant quantitative differences [3l[Tl]. Most importantly, 
this comparison highlights the fact that the HEA energy scale along the nodal direction 
is approximately twice that found in the hole-doped materials and the main valence 
band weight also lies at higher binding energy. 

To address the issues surrounding matrix elements that vary with incident photon 
energy, a 2D fitting procedure, described in [29], was used to extract the matrix 
elements along the nodal cut in both hole- and electron-doped materials. The 
momentum space profiles are shown in figure 3. In particular, the matrix elements 
significantly reduce the ARPES intensity upon approaching the F-point. However, 
taken together with the EDC and MDC dispersions for these materials, the matrix 
elements affect the appearance of the HEA, but can not account for the anomaly, also 
highlighted in a recent first-principles calculation [T8] . 



3. Theoretical Results 

The experimental findings clearly indicate the presence of a HEA in NCCO, not 
attributable to matrix element effects. This feature resembles the HEA found in 
both parent insulators [12] and hole-doped systems [3 [3l HI [5l [6l [7l [8] , and agrees well 
with the findings from previous studies on electron-doped compounds [Q] [101 [H] ■ This 
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section provides a unifying description of the underlying mechanism that gives rise to 
the HEA within the simple framework of the single-band Hubbard model. 

The two-dimensional, single-band Hubbard Hamiltonian serves as an effective, low 
energy model of the copper-oxide planes [301 F5T ] incorporating the effects of strong 
correlations. The single-band Hubbard Hamiltonian has the form 

ij,a i i 

where {i^ } is a set of tight-binding parameters where only nearest-neighbor t and next- 
nearest-neighbor t' are different from zero, c\^{ci^a) creates (annihilates) an electron 

with spin a at site i and rii^ — c\^Ci^a with Ui = Ui,-^ + Ui^i. The chemical potential 
H controls the electron filling; and, in what follows, the Hubbard repulsion U is set 
equal to the noninteracting electron bandwidth W = 8t. 

Here, the single-band Hubbard model is studied using the determinant quantum 
Monte Carlo technique [321 123]) an auxiliary- field method, to obtain the finite 
temperature, imaginary time propagator Gij (r) on a finite-size cluster with periodic 
boundary conditions. The details of the method can be found in [55] ■ The finite-size, 
square clusters used in this study have linear dimension N ^ 8 and the imaginary 
time interval, partitioned into L "slices" of size At = P/L, runs from to /3 = 1/T, 
the inverse temperature. The hopping t serves as the energy unit of the problem and 
unless otherwise noted At = f/f6t . 

The maximum entropy method |341 135] provides an effective means of obtaining 
the spectral function A(K,oj) from the imaginary time data on the grid {K} in 
momentum space. Once A(K,w) has been obtained, the single-particle self-energy 
I](K, uj) can be extracted using Dyson's equation and the tight-binding bandstructure. 
Assuming a weak momentum dependence to the self-energy, an interpolation routine 
provides the value of I](k, w) at an arbitrary point k in momentum space and Dyson's 
equation can be employed once again to compute ^(k, lj). 

The result from model calculations at half-filling provide a benchmark for the 
analysis and discussions of the behavior of the spectral function upon hole or electron 
doping. Figure 4 shows ^(k, w) for the single-band Hubbard model at half-filling 
(parameters are given in the figure caption). The canonical Hubbard bands are 
pronounced, but as in previous investigations [361 137] . low energy satellite bands 
near Ep also form. These branches of the lower and upper Hubbard bands serve 
as precursors for the quasi-particle bands that form upon hole or electron doping. 
These precursors also resemble the low energy features seen in ARPES experiments 
on the parent insulator [12] which provides a link for comparing the HEA in hole- and 
electron-doped systems that follows. 

Results from calculation of A(k, lj) for the single-band Hubbard model upon 
doping are shown in figures 5(a)-(d). The single-particle spectral function is shown for 
a representative system with parameters t' = —0.3 1 and /3 — 3/i at ~ 14% hole doping 
(/x = — 2.5i) in figures 5(a) and (b) while (c) and (d) show the results for a system 
with parameters t' = —0.2t and /3 = 3/i at ~ 16% electron doping (/x = 2.25 1). For 
each system the tight-binding parameters and even Fermi surface topology are quite 
similar. However, note the two sets of parameters are not simply related by a particle- 
hole transformation and, in general, the experimentally determined bandstructure 
parameters would not reveal such a symmetry between hole- and electron-doped 
materials. 
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Figure 4. Calculated single-particle spectral function A(k,Lu) at half-filling 
{fj, = 0.0 i), t' = — 0.3i and (3 = 3/t, along a nodal cut from (0,0) to (it, it). 
The inset shows a false-color plot of a portion of the same cut from (0, 0) to 
~ {tt/2, 7r/2) below the Fermi level, qualitatively similar to experimental results 
on the parent insulators |12) . The labels denote the lower Hubbard band (LHB), 
upper Hubbard band (UHB), and quasi-particle-like branches (QPB) that serve 
as precursors to the bands that form, and cross Ep, with hole or electron doping. 



For the hole-doped system the lower Hubbard band (LHB) is essentially localized 
near the F-point with a weak tail of spectral intensity extending toward the points 
(tt, tt) and (tt, 0). The decrease in spectral weight in these tails approximately coincides 
with the location in momentum space identified with the HEA. The weak appearance 
of the upper Hubbard band (UHB) results from hole-doping ^HIISS] as spectral weight 
is transferred into the LHB or more precisely the dispersive quasi-particle like band 
(QPB) that develops from low binding energy precursors found at half-filling [321 137] 
seen in figure 4. 

Along a nodal cut, the QPB is highly dispersive crossing the Fermi level at 
~ (7r/2,7r/2). Near (7r/4,7r/4) the spectral intensity in the band drops and the 
HEA appears as an apparent cross-over from the QPB to the LHB at an energy 
~ —0.5t to — 0.75i, the HEA energy scale. While the spectral intensity decreases 
as the QPB approaches the F-point, the coexistence of the LHB and QPB over a 
significant momentum interval supports the cross-over scenario presented here. 

In the electron-doped system, the dispersive QPB shown in figure 5(c) dips farther 
below the Fermi level than its hole-doped counterpart, to an energy ~ — 1.0 1 to — 1.5 1, 
in agreement with the comparison shown in figure 3. More importantly, the QPB arises 
from intermediate energy precursors in the UHB of the half-filled system [321 ISI] ; in 
contrast to the spectral weight forming the hole-doped QPB. While there appears to 
be a weak tail of the UHB at significantly higher energies (~ 3 1 to 4t), the dispersive 
feature and bulk of the UHB appear to meet in another "waterfall" well above the 
Fermi level. 

The appearance of the HEA as a "waterfall" at intermediate binding energies 
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Figure 5. Calculated single-particle spectral function A(k,a;) for ~ 14% hole- 
doping (fi = — 2.5t), t' = — 0.3t and /3 = 3/t in panels (a) and (b) and ~ 16% 
electron-doping (fi = 2.25 1), t' = — 0.2t and (3 = 3/t in panels (c) and (d). 
Panels (a) and (c) trace out the path (0, 0) — {tt,tt) — (tt, 0) — (0, 0) in the first 
BZ. Panels (b) and (d) focus on nodal cuts from ~ (—it/2, —n/2) to ~ (7'"/2, tt/2) 
where representative photoelectric matrix elements have been used in creating the 
intensity plots. The labels in panels (a) and (c) refer to the lower Hubbard band 
(LHB), upper Hubbard band (UHB), and quasi-particle band (QPB) with arrows 
provided to clarify the assignments. 



may be highlighted by applying photoelectric matrix elements, representative of 
experimentally derived values, to the results obtained from the calculations. Figures 
5(b) and (d) show A(k, w) multiplied by these matrix elements (see |2S]) along a 
nodal cut for hole- and electron-doped systems, respectively. In each case the matrix 
elements severely suppress the intensity near the F-point. In particular for hole- 
doped systems, this behavior is consistent with the interpretation of the QPB as a 
renormalized Zhang-Rice singlet (ZRS) band with photoelectric matrix elements that 
also would reflect the robustness of the Zhang- Rice picture in different regions of the 
BZ [30l|4^. This suppression of intensity, particularly in the hole-doped calculation, 
leads to the appearance of the HEA as a "waterfall" due to the overlap of spectral 
weight at intermediate energies for k away from the F-point. In the electron-doped 
calculation, the lack of significant spectral weight in the LHB leads to a simple 
suppression of intensity near the F-point without a significant overlap at intermediate 
energies. This result should be affected by the inclusion of additional valence band 
weight within multi-band treatments. 
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Figure 6. Schematic diagram for the multi-band structure of the copper-oxide 
piano highlighting the Zhang-Rice singlet (ZRS) band (blue) and upper Hubbard 
band (UHB) (red). The schematic also includes the multi-band lower Hubbard 
band (LHB), the Zhang-Rice triplet band (T) and the oxygen non-bonding band 
(NB). The ZRS band serves as the LHB in the single-band model, (a) Bands at 
half-filling highlighting the incoherent portions of the ZRS band (LHB for single- 
band) and UHB together with the quasi-particle band precursors. The chemical 
potential fi (Fermi level Ep) shifts as a function of (b) hole or (c) electron doping 
into the ZRS band (LHB for single-band) or the UHB, respectively. 



The relative spectral weight in the LHB and UHB as well as the HEA energy scale 
exemplify the dichotomy between hole- and electron-doped systems. A schematic of 
the band structure can be found in figure 6. In a multi-band treatment incorporating 
planar oxygen degrees of freedom, the UHB and ZRS band correspond roughly to the 
UHB and LHB in the single-band Hubbard model. Additional bands (proper multi- 
band LHB, Zhang-Rice triplet and non-bonding oxygen) serve as marker states and 
play no role in the single-band calculations. 

Upon hole doping in the single-band model, the chemical potential shifts into the 
LHB with a concomitant transfer of spectral weight from the UHB and increase in 
weight at energies near the chemical potential, forming the QPB (see figure 6(b) and 
changes to the multi-band ZRS band and UHB). In contrast, upon electron doping, 
the chemical potential shifts into the UHB with a transfer of weight to higher relative 
energy, again forming the bulk of the QPB (see figure 6(c)). The chemical potential 
and spectral weight shifts with doping are consistent with the treatment presented in 
Ref. [38]. Note with doping that the Mott gap does not simply collapse, but instead 
forms much of the intermediate energy regime below the HEA energy scale in the 
electron-doped calculation. This is in contrast to the hole-doped system in which the 
anomaly lies within the canonical LHB and the Mott gap has been effectively pushed 
above the Fermi level. 

The presented results support the conclusion that strong correlations and many 
body effects, here in the guise of the single-band Hubbard model, play a central 
role in the high energy anomaly. The calculations also echo some of the results of 
much earlier investigations into the spectral function of the single-band Hubbard 
model 133 ISZl mi- This study shows evidence for an anomaly in electron-doped 
materials and does not equate its energy scale, under either hole or electron doping, 
in any simple way with J that should be approximately equivalent in the two 
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types of materials. Band renornialization due to strong correlations incorporates 
coupling to paramagnon modes at these high energy scales, combined with all other 
renornialization pathways, to produce a transition, or cross-over, from a quasi-particle 
band at low binding energy near Ep to the incoherent lower Hubbard band, or more 
precisely oxygen valence bands, at higher binding energy. This behavior is consistent 
with experimental observations and not captured in weak coupling approaches that 
produce a simple "kink" in the dispersion at the HEA energy scale [Ml [TTl [18] . While 
similar conclusions have been reached in Ref. [10], the combined experimental and 
theoretical treatment presented here elucidates the complex interplay between doping, 
spectral weight transfer, and band renornialization that gives rise to the HEA beyond 
a straightforward chemical potential shift and a tight-binding analysis. 

Qualitatively, the doping and momentum dependence obtained from calculations 
are similar to the experimental results and, in many cases, minor adjustments to 
model parameters produce even surprisingly close quantitative agreement with the 
observed features [ll[31lll[51[ni[71[Hl[ni UHl [II]j especially those presented here on 
electron-doped NCCO. In the near nodal region, the agreement depends far less on 
the fine details of the underlying bandstructure; while moving toward the anti-nodal 
region, a quantitative agreement between theory and experiment depends heavily on 
the tight-binding fit parameters due to the already shallow band dispersion and relative 
strength of the renormalization, also demonstrated in the analysis of Ref. [10] ■ While 
a great deal of information has been obtained about the hole doping dependence of 
the HEA [21 [3] , the electron doping dependence remains largely unexplored to this 
point as do effects associated with changing matrix elements between different BZs. 
Preliminary results indicate that our theoretical simulations are qualitatively robust to 
changes in electron doping, as they are to hole doping, and also qualitatively reproduce 
dispersions in other BZs, but this remains the subject of ongoing investigation and 
detailed results will appear in a separate publication. It is also important to note that 
the single-band Hubbard model should only be viewed as a low energy effective model 
of the cuprates. A better understanding of the role played by strong correlations would 
come from careful analysis of calculations explicitly incorporating separate copper and 
oxygen degrees of freedom that would capture the formation of Zhang- Rice singlets 
[30l [40] and the spectral weight associated with the valence bands, assumed to have 
strong oxygen character. 
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